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Resumen 


MORALES-VELÁZQUEZ, M.I., APARICIO, J. ée VALDÉS, 
J.B. Pronóstico de avenidas utilizando el filtro de Kalman 
discreto. Tecnología y Ciencias del Agua. Vol. V, núm. 2, marzo- 
abril de 2014, pp. 85-110. 


Se evalúa la utilidad y aplicabilidad del algoritmo del filtro de 
Kalman discreto en la predicción de caudales a corto plazo. 
El algoritmo se aplica a la cuenca propia de la presa Ángel 
Albino Corzo (Peñitas), parte del Sistema Hidroeléctrico 
Grijalva, y a la estación hidrométrica Sayula. El algoritmo 
se utiliza para determinar la función de respuesta en la 
cuenca y con ello pronosticar los caudales de entrada al 
embalse. Para esto se usan tanto los registros de caudal y de 
precipitación provenientes de las estaciones climatológicas 
ubicadas en la cuenca como los caudales calculados con el 
tránsito inverso en el vaso. En el análisis se evalúan varios 
intervalos de tiempo de pronóstico, así como diferentes 
tipos de función de respuesta y parámetros asociados. Los 
resultados son evaluados por medio del coeficiente de Nash- 
Sutcliffe, obteniéndose valores muy aceptables, de manera 
que el filtro se considera aplicable al pronóstico de avenidas 
a corto plazo, destacando su utilidad como una herramienta 
de apoyo en el desarrollo de políticas de operación y control 
de los embalses. 


Palabras clave: pronóstico de avenidas a corto plazo, filtro 
de Kalman, Peñitas, recursivo. 


Introducción 


A lo largo de la historia y en particular en 
años recientes, se han presentado importantes 
inundaciones en México y otros lugares en el 
mundo. Por ejemplo, en el sureste mexicano 
destacan las ocurridas en Chiapas en 1998; en 
Quintana Roo, Chiapas y Veracruz en 2005, 
como consecuencia de los huracanes Emily, Stan 
y Wilma; y las de Tabasco en 2003, 2010 y 2007 
(Salas y Jiménez, 2004). Esta última, en particular, 
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provocó el aumento de los niveles del agua 
en el sistema de ríos de la cuenca del Grijalva, 
causando inundaciones muy significativas en la 
mayor parte del estado de Tabasco y en Chiapas 
(Aparicio et al., 2009; Rivera-Trejo et al., 2009). 
Para hacer frente a estos episodios de inun- 
daciones mediante sistemas de alerta y criterios 
de operación de la infraestructura de control 
es importante contar con sistemas confiables y 
precisos de pronóstico de avenidas. El propósito 
del presente trabajo es evaluar la utilidad del 
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algoritmo del filtro de Kalman discreto para el 
pronóstico de avenidas a corto plazo, tomando 
como caso de estudio la cuenca propia de la 
presa Ángel Albino Corzo (Peñitas) (figura 1). 


Metodología 


Se utilizará el filtro de Kalman (Kalman, 
1960) para identificar la función de respuesta 
de la cuenca y con ella los gastos de entrada 
al vaso de la C.H. Peñitas mediante la con- 
volución de dicha función de respuesta con 
las precipitaciones, o con los gastos y preci- 
pitaciones registrados en un lapso anterior al 
tiempo analizado. 

El filtro de Kalman (Kalman, 1960) describe 
una solución recursiva para el problema del 
filtrado lineal de datos discretos. El filtro 
trata de estimar el estado xE%R"” (espacio 
n-dimensional de los números reales) en el 
tiempo kAt (en adelante k) de un proceso 
controlado en tiempo discreto, descrito por la 
ecuación matricial diferencial lineal estocástica 
(Welch y Bishop, 2001): 


Xy = AXz 1 +BU¿ +01 (1) 


con una medición z € R”, que está representada 
por: 


2, = Hx, +0, (2) 


La matriz A, en la ecuación (1) relaciona 
el estado en el tiempo previo k-1 al estado que 
ocurre en el momento k; esta matriz puede 
cambiar en el tiempo, pero habitualmente se 
considera constante (Welch y Bishop, 2001). 
La matriz B, , relaciona el control opcional de 
entradas u E %' con el estado x y la matriz H 
en la ecuación (2) relaciona el estado con la 


mxn 


medición z,. En la práctica en hidrología, B 


nxl 


por lo general también se considera constante 
(Ibid.). 

Por otra parte, las variables aleatorias 10, 
y v, en las ecuaciones (1) y (2) representan el 
ruido en el proceso y el ruido en la medición, 
respectivamente. Se acepta que estas variables 
son independientes una de otra, que son ruido 
blanco (es decir, sus valores de señal en dos 
tiempos diferentes no guardan correlación 
estadística) y que tienen una función de 
distribución de probabilidad normal (Drécourt, 
2003; Kottegoda, 1980; Kim et al., 2004): 


p(w)-N (0,5) (8) 


p(v)-N(0,R) (4) 


dl 


Vaso superior 


Presa Malpaso 723.6 km? 


Presa Peñitas 


Vaso inferior 
532.9 km? 


e Estaciones 
hiddroclimatológicas 


Figura 1. Cuenca de la C.H. Peñitas (Aguilar et al., 2009). 
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S representa la matriz de covarianza de la 
perturbación del proceso y R la matriz de covarianza 
de la perturbación de la medición. Ambas matrices 
podrían cambiar en el tiempo; sin embargo, por 
simplicidad, en la práctica también se suelen 
suponer constantes y pueden ser representadas 
a su vez de la siguiente manera (Simon, 2001): 


S= Etoo; | (5) 
R= El ojo; | (6) 


donde El+] es la esperanza matemática. Estas 
matrices representan, respectivamente, la 
variación de los errores o perturbaciones en las 
estimaciones del proceso o sistema analizado, y 
de las mediciones del mismo con respecto a la 
media del proceso. 

El filtro de Kalman utiliza un control de 
retroalimentación, ya que estima el proceso en 
un momento en el tiempo y entonces obtiene la 
retroalimentación por medio de la medición de 
los datos observados (Welch y Bishop, 2001). 

Desde este punto de vista, las ecuaciones 
que se utilizan para derivar el filtro de Kalman 
se dividen en dos grupos (Ibid.): 


e Ecuaciones de actualización en el tiempo o 
ecuaciones de predicción. 

e Ecuaciones de actualización con los datos 
observados. 


Las del primer grupo son responsables 
de proyectar hacia adelante en el tiempo el 
estado actual, tomando como referencia el 
estado en el momento previo k — 1, así como 
de actualizar las estimaciones de la covarianza 
del error, para obtener las estimaciones a priori 
del estado para el próximo paso en el tiempo. 
El segundo grupo de ecuaciones es responsable 
de la retroalimentación, es decir, de incorporar 
nueva información dentro de la estimación 
anterior, con lo cual se llega a una estimación 
mejorada del estado a posteriori. Así, el filtro de 
Kalman constituye un algoritmo de proyección- 
corrección al pronosticar el nuevo estado y su 
incertidumbre, y corregir la proyección con la 
nueva medición. 


La aplicación del filtro de Kalman podría 
resumirse en los siguientes pasos (Welch y 
Bishop, 2001), mostrados en la figura 2; el acento 
circunflejo indica estimación y el superíndice 
denota que la estimación es a priori: 


1. Generar un pronóstico del estado hacia 
adelante en el tiempo, tomando en cuenta 
toda la información disponible hasta ese 
momento, partiendo de la proposición 
inicial del estado y de la matriz de 
covarianza del error Pk. 

2. Corregir la proyección del estado (proceso 
de actualización con las mediciones). 
La primera tarea dentro de esta etapa 
es el cálculo de la ganancia de Kalman, K.. 
Este factor de ponderación o ganancia se 
selecciona de tal forma que minimice la 
covarianza del error de la nueva estimación 
del estado, teniendo en cuenta tanto el 
error probable de las mediciones como la 
incertidumbre sobre la representación del 
estado del sistema. 

3. Medir el proceso para obtener z, y entonces 
generar una nueva estimación del estado, 
incorporando la nueva medición. 

4. Obtener una nueva estimación de la matriz 
de covarianza del error, para después 
valorar de nuevo el estado. 


Después de cada par de actualizaciones, 
tanto del estado como de la medición, el 
proceso se repite, partiendo de las nuevas 
estimaciones del estado y de la covarianza del 
error. Esta naturaleza recursiva es una de las 
características más importantes del filtro de 
Kalman. 


Aplicación del filtro de Kalman discreto 
(FKD) a la cuenca propia de la C.H. 
Peñitas 


Para la predicción de los caudales es necesario 
contar con hietogramas de precipitación efec- 
tiva. Para esto se obtienen las curvas que relacio- 
nan el índice de precipitación antecedente (IPA) 
y el índice de infiltración media (p) (Aparicio, 
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CS => 


Actualización en el tiempo 
(“Predicción”) 


(1) Proyección del estado hacia adelante 


X = AX, +Bu, 


+ Proyección de la covarianza del error 
hacia adelante 


Po = AP, ¡+47 +S 


W 


Estimaciones iniciales para X;1 
y para Py 


Actualización con la medición 
(“Corrección”) 


(2) Cálculo de la ganancia de Kalman 
T ya cl 
K, = PH" (HP¿H" +R) 


(3) Actualización de la estimación mediante 
la nueva medición 


Xi. =* +K (2, -Hí;) 
(4) Actualización de la covarianza del error 


P¿ (1-K,H) Pz 


EX 


Figura 2. Diagrama completo de la operación del filtro de Kalman (Welch y Bishop, 2001, p. 24). 


2011). Esta relación permite deducir, a partir 
de la lámina de lluvia real, la lluvia efectiva 
que actuará como estímulo en la cuenca, 
produciendo los escurrimientos de entrada. 

En este trabajo, el análisis se realiza consi- 
derando los registros de precipitación y de 
escurrimiento en las estaciones ubicadas en la 
cuenca inferior (cuadro 1). El periodo analizado 
es del 31 de octubre de 2005 a las 24:00 h al 13 
de febrero de 2012 a las 24:00 horas. 

Para la construcción de la relación IPA vs. 
q, la precipitación es la media hp obtenida 
mediante polígonos de Thiessen, mientras que 
los gastos son los registrados en la estación 
hidrométrica Sayula. La curva IPA vs. q 
resultante es (ver figura 3): 


q =-4.213 In(IPA) +24.877 (7) 


Con la ecuación (7) se calculó la lluvia 


efectiva hp de la siguiente manera: 


efectiva 


Si p> hp > WPetectiva =0 (8) 


Si p< hp q HPetectiva 5 hp-p (9) 

No existen suficientes estaciones hidromé- 
tricas disponibles en la cuenca de la C.H. 
Peñitas para determinar el caudal real de 
entrada al embalse. Como se verá más adelante, 
el filtro se evaluó utilizando los caudales 
medidos en la estación hidrométrica Sayula 
y, para el pronóstico de caudales en la C.H. 
Peñitas, se utiliza como alternativa el tránsito 
inverso modificado en vasos o antitránsito, 
mediante el cual, conociendo las salidas totales 
en la cuenca, así como la evolución de los 
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Cuadro 1. Estaciones climatológicas e hidrométricas utilizadas. 


Climatológicas Hidrométricas 
Cc Rómul Juan del 
RO A Zapata | Tzimbac NO Ocotepec Malpaso Grijalva | Aza-Pac | Tzimbac | Aza-Pac 
superior Calzada Eupeñor 
Cuenca P0n E > 
, , Peñitas[Ocotepec| Zapata | Sayula | Juan del Grijalva Superior Sayula 
inferior 
respuesta. Para ejemplificar el proceso de 
IPA vs y cálculo, supóngase que se tiene un hietograma 
Subcuenca inferior. Río Sayula e $ . 
12 > de precipitación efectiva como el mostrado en 
Ss . . 
- 1 a y =-4.213 Ins) + 24.877 la figura 5a, el cual produce un hidrograma 
Ñ R?= 0.7433 e y 
3 8 == como el mostrado en la figura 5c a través del 
Es * hidrograma unitario instantáneo presentado 
E + e ] 
+ 
24 + == en la figura 5b. 
3 á . Durante el pronóstico de una avenida, las 
£ ] barras del hietograma no son conocidas en su 
0 50 100 150 200 250 


Índice de precipitación antecedente IPA (mm) 


Figura 3. Curva IPA us q. 


niveles en el vaso, se estiman los volúmenes 
o gastos totales de entrada al mismo (Aldama 
y Aguilar, 1996). Los cálculos del antitránsito 
para la cuenca inferior se considerarán, dentro 
del contexto del FKD, como las mediciones de 
los ingresos a la C.H. Peñitas para los efectos de 
este trabajo, ya que contienen implícitamente el 
comportamiento de toda la cuenca (figura 4). 

El cálculo del antitránsito presenta oscila- 
ciones significativas en las entradas, situación 
común en este tipo de análisis (Aldama y 
Aguilar, 1996). Para lograr una señal con 
menos ruido, se aplican promedios móviles de 
los cinco últimos valores medidos (Aguilar et 
al., 2009). 


Algoritmo del FKD utilizando el hidrograma 
unitario instantáneo como función de 
respuesta 


En esta modalidad, el cálculo de los caudales 
se hizo mediante la aplicación del hidrograma 
unitario instantáneo (HUI) como función de 


totalidad, pero se sabe la función de respuesta; 
es decir, el HUL, cuyas ordenadas se pueden 
representar como un incremento entre el gasto 
en un instante y el siguiente (como en el caso de 
3q,, en la figura 5b); los gastos del hidrograma 
de escurrimiento directo estarán representados 
como: 


Q; =hp9y 
Q) =Q, + hp,9q,> + hp,qu 


Qs =05 +/p,0955 + hp25954 + hp30q 43 + 1p4532 
+/1p50,1 + hp¿Sq, 


Si se desea predecir cada uno de los 
caudales con base en lo registrado en k — 1, los 
valores resaltados en negritas en el sistema 
anterior representan un error en la predicción 
de los mismos, debido a que aún no se registra 
la precipitación correspondiente al tiempo k. 

Con base en lo anterior, y de forma gene- 
ralizada, los caudales pueden representarse 
como: 


Q, =Hx,+0Q, ¡+8 (10) 
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Descargas vertedor Peñitas 
C.H. Peñitas “AR 


Descargas canal JDG + Túnel 
ade Caído JDG ) 
€ | mel 
Pa eo y ECP (Peñitas) 
A *. | Infitración 
ECP (DG) ECP (JDG) 7 Cuenca inferior (Peñitas) 
¡ ( io Río principal Sayula ) 
( Infiltración Cuenca superior (Juan del Grijalva) E ls ESE 
AAN y Río principal Tzimbac PERA SAA P o 
ps 7 <e > 1 A: »] 
DoS o » aporación SGL 5 / 


e 


Valpa 297 


b) 


Figura 5. Hidrograma unitario instantáneo. Ejemplo. 


Á 

o 

% Q, = caudal registrado en el momento f. 

mita x, = vector que contiene los incrementos entre 

eS n-i-1,n-i 7 

a . 

o ordenadas sucesivas del HUI. 

E H = vector con las precipitaciones anteriores 

a H=|hp, 1,Mp,_2,..Mp,_n ]; x,= 974-211 al tiempo analizado, disponibles en el 

E 6 hietograma de precipitación efectiva. 

> E E A 

$ 4-1, n = cantidad de barras en el hietograma de 

e Seca ., . 

S | 9 m0+1 precipitación efectiva. 

=> e =error al predecir el caudal, absorbido 
donde: por la parte en la convolución en la que 
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intervienen las precipitaciones que aún no 
son registradas. 


La variable a estimar mediante la aplica- 
ción del algoritmo del FKD con base en lo 
desarrollado, serán los incrementos entre 
las ordenadas de la función de respuesta de 
la cuenca; en este caso, representada por el 
HUI. Estos incrementos son entendidos en la 
nomenclatura del filtro por el término estado, 
denotado por *,. 

Debido a que no se conoce con exactitud 
el valor que tomará el estado inicial, su valor 
esperado se considera de manera arbitraria 
como un vector nulo [0], , de dimensión n, 
donde n representa el número de barras del 
hietograma de precipitación efectiva ante- 
riores al tiempo analizado, consideradas en 
la convolución con los incrementos entre las 
ordenadas de la función de respuesta para 
obtener los caudales de entrada: 


La = [0]... (11) 
mientras que la matriz de covarianza del error 
del estado inicial será: 

Pr =09n! (12) 
donde n es un escalar lo suficientemente 
grande que refleja la incertidumbre de los 
valores supuestos para el estado inicial, en 


concordancia con Valdés et al. (1980); I es la 
matriz identidad de dimensiones n x n. 


Ecuaciones de pronóstico del estado 


Proyección del estado hacia adelante 
(pronóstico): 


%¿ = AX, ¡ +BUj».. (13) 


A, es una matriz identidad que relaciona el 
estado en el momento k-— 1 con el que ocurre en 
k; Ne es el estado estimado a priori; B,_ estará 
definida como una matriz identidad con las 


mismas dimensiones que la matriz A; mientras 


que u,, que representa el control opcional de 
entradas, será considerado como un vector 
de ceros de dimensiones n x 1 debido a que el 
sistema no contiene variables que se puedan 
controlar y que influyan en su respuesta. 


Proyección de la matriz de covarianza del 
error hacia adelante (P;): 
Po = AP, ¡A +S... (14) 

S es la matriz de covarianza de la pertur- 
bación del proceso de dimensiones n x n. En 
este apartado se considera de tres formas dis- 
tintas (como una matriz de ceros, una matriz 
cuya diagonal principal depende de los errores 
entre el estado pronosticado y actualizado, o 
bien como un valor muy pequeño cercano a 
cero), como se verá más adelante. 

Con las ecuaciones anteriores es posible 
obtener el primer pronóstico del estado, así 
como de la matriz de covarianza del error. 
Una vez obtenida la función de respuesta, el 
pronóstico de los caudales de entrada al vaso de la 
C.H. Peñíitas se calcula como sigue: 

Qest1 + H-Xz +Q61 (15) 

Donde X; contiene los incrementos entre 
las ordenadas (q) del HUL obtenidas con 
el FKD; Q,,, es el caudal pronosticado en 
la cuenca propia de la C.H. Peñitas; Q, , es 
el caudal registrado en el tiempo k — 1, y 
Has =[ Hp, 1, Hpy-21=»>+ HP] es un vector 
que contiene los valores a convolucionar con 
la función de respuesta, conformado por las 
mediciones de precipitación efectiva previas al 
momento k analizado. 


Actualización o corrección del pronóstico 
Cálculo de la ganancia de Kalman (K,) 

Para el caso del presente trabajo, la matriz de 
covarianza de la perturbación de la medición 


¡y donde 
a. es una constante de proporcionalidad que 


R (paso 2 de la figura 2) será R=a. + Q 
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representa un error constante igual a una 
fracción del caudal medido en el tiempo anterior 
k-1 (Valdés et al., 1980). En la actualización del 
pronóstico con base en las mediciones actuales 
(paso 3 de la figura 2), z, será, con base en la 
ecuación (10), el caudal real medido en kmenos 
el caudal medido en k- 1, tomado directamente 
del cálculo del antitránsito o del registro de la 
estación hidrométrica, según sea el caso. 

Una vez llevada a cabo la actualización 
del pronóstico de la función de respuesta en 
la cuenca, se actualiza el pronóstico de los 
caudales de entrada al vaso de la C.H. Peñitas 
de la siguiente manera: 

Qest2 =H La +Qí1 (16) 
donde %, son los incrementos entre las orde- 
nadas de la función de respuesta corregidos o 
bien actualizados. 

Dentro de la rutina de aplicación del FKD, 
se obtuvieron los siguientes estadísticos que 
permiten medir de forma cuantitativa la pre- 
cisión del pronóstico: 


+ Coeficiente de Nash-Sutcliffe (Nash y 
Sutcliffe, 1970). 

e Media y desviación estándar tanto para la 
serie de caudales originales como para los 
caudales pronosticados y actualizados. 

e  Coeficientes de correlación entre la serie 
original de los caudales de entrada a la 
C.H. Peñitas, calculados con el antitránsito 
y la serie de caudales pronosticados y 
actualizados con el FKD. 


Pronóstico de caudales en la cuenca completa 
de la C.H. Peñíitas, utilizando el HUI como 
función de respuesta 


Con fines experimentales, se evaluaron tres 
diferentes formas de cálculo de la matriz de 
covarianza de la perturbación o ruido en el 
proceso S. En cada una de ellas se probaron 
diferentes valores de 11, a. y N, pero por razones 
de espacio sólo se presentan los resultados 


obtenidos con la mejor combinación, usando 
un intervalo de pronóstico de 1 h con fines de 
evaluación. 


Aplicación del FKD usando S =[0],. 


Para este caso, la matriz S se obtuvo consi- 
derando que está conformada por una matriz 
nula de dimensiones n x n. Los mejores resul- 
tados fueron obtenidos considerando 3 h 
previas de precipitación efectiva en la matriz 
H. En la figura 6a se muestra la precipitación 
efectiva media calculada y en la figura 6b se 
tienen los caudales observados (calculados con 
antitránsito) y pronosticados (FKD). La línea 
horizontal entre las 1.026 x 10* y 1.796 x 10* h 
indica ausencia de registros de caudal. En la 
figura 7 se presentan los gastos pronosticados 
contra los observados. La línea diagonal a 45” re- 
presenta una coincidencia perfecta. La relación 
entre gastos observados contra actualizados 
(no presentado en este trabajo) resulta similar. 
En las figuras 8, 9 y 10 se presentan algunas 
avenidas específicas dentro del lapso analizado 
para una comparación detallada. Se usó un 
solo valor de gasto anterior para el pronóstico, 
observándose una coincidencia notable en los 
hidrogramas mencionados, en particular en los 
gastos máximos. 

En el cuadro 2 se resumen los estadísticos 
utilizados para evaluar cuantitativamente la 
precisión de los pronósticos tomando a = 0.3, 
N=1000, n =3. 

La media y desviación estándar de los cau- 
dales observados fueron 759.26 m*/s y 428.79, 
respectivamente, mientras que esos valores 
para los gastos pronosticados fueron de 760.04 
y 429.58 m*/s. Se observa que se preservan de 
forma razonable estos estadísticos. 


Aplicación del FKD considerando S = 0 


Con el objeto de realizar un análisis de 
sensibilidad, se probaron valores diferentes de 
cero, pero pequeños para las componentes de la 
diagonal principal de la matriz S, manteniendo 
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Registros de precipitación efectiva horaria 
Cuenca propia de la C. E Peñitas (2005- a 


a TI 
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Figura 6. Registros de precipitación y caudal utilizados en el análisis de la cuenca propia de la C.H. Peñitas 


y pronóstico de caudales con 1 h de anticipación. 


Caudales “observados” vs. caudales pronosticados 


Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 7. Caudales observados versus pronosticados en la cuenca propia de la C.H. Peñitas, considerando el HUI 


como función de respuesta y S =[0]. 


la ecuación (3). La respuesta fue evaluada 
considerando los mismos valores de a, N y 
n, que en el análisis anterior. Los resultados 
se muestran en el cuadro 3. Se observa una 
mejoría marginal entre los resultados. Por lo 
anterior, no se recomienda usar esta opción. 


Aplicación del FKD considerando S =X,-X, 


Una alternativa adicional fue obtener el valor 
de S como una diferencia entre las ordenadas 
de la función de respuesta pronosticada X, y 
actualizada X,, como se describe en Kim et al., 
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Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 8. Caudal observado y pronosticado con el FKD (avenida A del 23/11/2007 01:00 h al 28/11/2007 01:00 h), 
aplicado a la cuenca propia de la C.H. Peñitas, considerando el HUI como función de respuesta y S = [0]. 
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Figura 9. Caudal observado y pronosticado con FKD (Avenida B del 25/09/2008 03:00 h al 26/09/2008 19:00 h), 
aplicado a la cuenca propia de la C.H. Peñitas, considerando el HUI como función de respuesta y S = [0]. 
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Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 10. Caudal observado y pronosticado con FKD (Avenida C del 31/10/2009 9:00 h al 03/11 /2009 15:00 h), aplicado a la 
cuenca propia de la C.H. Peñitas, considerando el HUI como función de respuesta y S = [0]. 


Cuadro 2. Resumen de estadísticos para la aplicación del FKD, considerando como función de respuesta el HUI 


y S=10] 
Concepto Valor 
Nash-Sutcliffe (Q “observado” - Q pronosticado) 0.977353 
Nash-Sutcliffe (Q “observado” - Q actualizado) 0.977534 
Correlación (Q “observado” - Q pronosticado) 0.99727 
Correlación (Q “observado” - Q actualizado) 0.99729 


Cuadro 3. Resumen de estadísticos para la aplicación del FKD, considerando como función de respuesta el HUI y S =0. 


A Nash-Sutcliffe Nash-Sutcliffe 
(Q “observado” — Q pronosticado) (Q “observado” - Q actualizado) 
0.1 0.977031 0.978011 
0.001 0.977389 0.977709 e 
0.0001 0.977412 0.977633 y 
0.00001 0.977389 0.977581 z 
0.0000001 0.977354 0.977535 S 
5 
2004. Entonces la matriz S estará compuesta ticada y la actualizada, acumuladas en un > 
en su diagonal principal por la desviación de lapso anterior Af. Los mejores resultados 5 
las diferencias entre los incrementos de las fueron obtenidos considerando 3 h previas El 
ordenadas de la función de respuesta pronos- de precipitación efectiva en la matriz H y un e 
3 
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lapso acumulado de los errores en la función 
de respuesta de 16 pronósticos previos. En la 
figura 11 se muestran los caudales observados 
contra los pronosticados. Se observan, en este 
caso, algunas diferencias importantes, así 
como también en la figura 12, donde se mani- 
fiesta una mayor dispersión respecto a los 
casos anteriores. Estas diferencias se aprecian 
también en las avenidas aisladas de las figuras 


13, 14 y 15. Por consecuencia, los valores del 
coeficiente de Nash-Sutcliffe y los coeficientes 
de correlación son menores, como se muestra 
en el cuadro 4. 

La media y desviación estándar de los 
caudales observados fueron de 759.26 y 428.79 
m?/s, respectivamente, mientras que esos valo- 
res para los gastos pronosticados fueron de 
760.35 y 432.24 m/s. 


Caudales “observados” y pronosticados 


eb Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 11. Caudales observados y pronosticados en la cuenca de la C.H Peñitas, considerando el HUI 


como función de respuesta y S=X,—X,. 


Caudales “observados” vs. caudales pronosticados 


Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 12. Caudales observados versus pronosticados en la cuenca de la C.H Peñitas, considerando el HUI 


como función de respuesta y S=X,—X,. 


Morales-Velázquez et al., Pronóstico de avenidas utilizando el filtro de Kalman discreto 


500 


400 


300 


Caudal Q (mé /s) 


200 


100 


PA A A AAA 


Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 


—+— Caudales “observados” (calculados con antitránsito) 
——-— Caudales pronosticados (FKD) 


=E 


de 


1.805 1.8062 1.8074 1.8086 1.8098 1.811 1.8122 1.8134 1.8146 1.8158 


1.817 
Tiempo (horas) x 10% 


Figura 13. Caudal observado y pronosticado con FKD (avenida A del 23/11 /2007 01:00 h al 28/11/2007 01:00 h), 


aplicado a la cuenca propia de la C.H. Peñitas, considerando el HUI como función de respuesta y S =X,-X.. 


2 000 


1800 


1600 


1 400 


1200 


1 000 


Caudal Q (m?/s) 


Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 


—4— Caudales “observados” (calculados con antitránsito) : 
——-— Caudales pronosticados (FKD) A 


Lp FITTTOTTITITITITTTATTNTTTT Ae A A ana al 


¡ | | ¡ | 
0 

2.542 2.5426 2.5432 2.5438 2.5444 2.545 2.5456 2.546 

Tiempo (horas) x 10% 


Figura 14. Caudal observado y pronosticado con FKD (avenida B del 25/09/2008 03:00 h al 26/09/2008 19:00 h), 


aplicado a la cuenca propia de la C.H. Peñitas, considerando el HUI como función de respuesta y S =X —X'. 
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Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 15. Caudal observado y pronosticado con FKD (avenida C del 31/10/2009 09:00 h al 03/11/2009 15:00 h), 
aplicado a la cuenca propia de la C.H. Peñitas, considerando el HUI como función de respuesta y S =X,-X.. 


Cuadro 4. Resumen de estadísticos aplicando el FKD a la cuenca propia de la C.H. Peñitas, considerando como función de 


respuesta el HUI y S=X,-X, acumulado de 16 pronósticos previos. 


Concepto Valor 
Nash-Sutcliffe (Q “observado” - Q pronosticado) 0.977298 
Nash-Sutcliffe (Q “observado” - Q actualizado) 0.977709 
Correlación (Q “observado” - Q pronosticado) 0.99726 
Correlación (Q “observado” - Q actualizado) 0.99731 


Discusión de resultados Algoritmo del FKD usando una función 
de respuesta modificada 

Al evaluar las tres diferentes alternativas de 
cálculo de la matriz S, se observa que cuando 


la matriz S es nula, los valores del coeficiente a eri 
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de Nash-Sutcliffe son muy aceptables (NS = 
0.97735), lo que significa que el pronóstico se 
realiza con buena aproximación. Al usar la 
matriz S, compuesta por números cercanos a 
cero, el valor del coeficiente de Nash-Sutcliffe 
(NS = 0.97741) mejora sólo de modo marginal. 
Si dicha matriz se conforma con los errores 
entre la función de respuesta pronosticada 
y actualizada, los resultados son menos 
satisfactorios. Por lo tanto, la mejor alternativa 
de cálculo es cuando S = [0], por lo que ésta 
será utilizada en lo subsecuente. 


que la función de respuesta predicha estaría 
representada por el hidrograma unitario instan- 
táneo y se usó exclusivamente la precipitación 
efectiva en la convolución con la función de 
respuesta. En lo que sigue, se evalúa una 
segunda alternativa, que incluye, además de 
la precipitación efectiva, los caudales medidos 
en un lapso anterior al tiempo analizado en la 
matriz H. El periodo de análisis considerado es 
el mismo que en el pronóstico anterior, es decir, 
del 31 de octubre de 2005 a las 24:00 h al 13 de 
febrero de 2012 a las 24:00 h. 
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Las variables incluidas en la aplicación del 
algoritmo del FKD representan básicamente lo 
descrito en el apartado anterior, salvo que en 
este caso las dimensiones de las matrices A, H, 
S, etc., son n + nQ, en lugar de n donde n y nQ 
son el número de episodios de precipitación 
efectiva y de caudales, respectivamente, 
considerados en la convolución realizada con 
la función de respuesta. Una vez obtenida la 
función de respuesta con la aplicación del FKD, 
el pronóstico de los caudales de entrada al vaso de 
la C.H. Peñitas se hace a través de la ecuación: 

Qos = HXz (17) 
donde X; son las ordenadas de la función de 
respuesta pronosticada y H,,, y, =[Q,y Qyy 
Quer Pro HP yor + Hp, ,] será un vector que 
contiene los valores de caudal 10 y precipita- 
ción efectiva n utilizadas en la convolución 
con la función de respuesta pronosticada para 
obtener los caudales de entrada al vaso. 

Una vez llevada a cabo la actualización del 
pronóstico de la función de respuesta en la 
cuenca (ecuaciones de actualización), siguien- 
do el proceso descrito con anterioridad, se 
actualiza el pronóstico de los caudales de 
entrada al vaso de la C.H. Peñitas en la forma: 

Queso =H Xy (18) 
donde Y, es la función de respuesta corregida o 
bien actualizada. 


Pronóstico de caudales en la estación 
hidrométrica Sayula usando la función de 
respuesta modificada 


Con el objeto de evaluar la aplicabilidad del 
FKD con la función de respuesta modificada, 
se usaron los registros de caudal en la estación 
hidrométrica Sayula (figura 1), considerando 
un periodo de análisis del 01/01/2010 al 
14/02/2012. 

Los mejores resultados se obtuvieron con- 
siderando 48 h previas de lluvia efectiva y 
1 h de medición previa de caudal, con un 


intervalo de pronóstico de 1 h. En la figura 
ló6a se muestra la precipitación efectiva media 
registrada en la estación hidrométrica Sayula, 
mientras que en la figura 16b se presentan los 
caudales observados y pronosticados mediante 
la aplicación del FKD. En la figura 17 se 
presentan los gastos pronosticados contra los 
“observados”. La línea a 45” representa una 
coincidencia perfecta, pudiéndose notar que en 
algunos de los casos el caudal pronosticado es 
subestimado y en otros sobreestimado, pero en 
general los resultados son muy satisfactorios. 
La relación entre gastos observados contra 
actualizados no presentados aquí resulta 
similarmente satisfactoria. Las figuras 18, 19, 
20 y 21 presentan algunas avenidas aisladas 
dentro del periodo de estudio para un mejor 
análisis. 

En el cuadro 5 se presentan los estadísticos 
utilizados para evaluar cuantitativamente la 
precisión de esta alternativa de aplicación del 
FKD. Los valores de a = 0.3 y N = 1 000 son 
iguales a los utilizados en los análisis anteriores 
y el valor de la matriz S = O corresponde al 
mejor resultado obtenido de la evaluación de 
las alternativas de cálculo de $. 

La media y desviación estándar de los 
caudales observados fueron de 21.51 y 24.37 
m/s, respectivamente, mientras que esos 
valores para los gastos pronosticados fueron 
de 21.53 y 24.61 m*/s. De nuevo se preservan 
razonablemente bien estos estadísticos. 


Discusión de resultados 


Los resultados obtenidos al aplicar el algoritmo 
del FKD a la estación hidrométrica Sayula son 
muy aceptables. Se encontró que el intervalo 
de 48 h previas de precipitación representa 
la mejor combinación en el pronóstico de los 
caudales. 


Pronóstico de caudales en la C.H. Peñitas 
usando la función de respuesta modificada 


Se presentan los resultados de la aplicación del 
filtro con la función de respuesta modificada 
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a) Registros de precipitación efectiva horaria 
Estación hidrométrica Sayula. Análisis horario (2010-2012) 
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Figura 16. Registros de precipitación y caudal utilizados en el análisis para la estación hidrométrica Sayula, incluyendo 
pronóstico de caudales, considerando la función de respuesta modificada. 


Caudales observados vs. caudales pronosticados 
Estación hidrométrica Sayula. Análisis horario (2010-2012) 
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Figura 17. Caudales observados versus caudales pronosticados de la estación hidrométrica Sayula, considerando 


la función de respuesta modificada. 
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usando como medición los caudales calculados Análisis horario 

con antitránsito para intervalos de pronóstico 

de 1 y 24 h. El lapso utilizado fue del 31 de Los mejores resultados se logran conformando 
octubre de 2005 al 13 de febrero de 2012. la matriz H por las mediciones de 12 h previas 
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Caudales observados y oronosticados 
Estación hidrométrica Sayula. Análisis horario (2010-2012) 
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Figura 18. Caudal observado y pronosticado con FKD (avenida A del 24/06/2010 al 30/06/2010) para la estación hidrométrica 


Sayula, considerando la función de respuesta modificada. 


Caudales observados y pronosticados 
Estación hidrométrica Sayula. Análisis horario (2010-2012) 
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Figura 19. Caudal observado y pronosticado con FKD (avenida B del 3/07/2010 al 13/07/2010 coincide con el huracán Alex) 
para la estación hidrométrica Sayula, considerando la función de respuesta modificada. 


Tecnología y 


Ciencias del Agua, vol. V, núm. 2, marzo-abril de 2014 


Morales-Velázquez et al., Pronóstico de avenidas utilizando el filtro de Kalman discreto 


Caudales observados y pronosticados 
Estación hidrométrica Sayula. Análisis horario (2010-2012) 
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Figura 20. Caudal observado y pronosticado con FKD (avenida C del 16/09/2010 al 12/09/2010 coincide con el huracán Karl) 
para la estación hidrométrica Sayula, considerando la función de respuesta modificada. 
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Figura 21. Caudal observado y pronosticado con FKD (avenida D del 25/09/2010 al 02/10/2010) para la estación hidrométrica 
Sayula, considerando la función de respuesta modificada. 
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Cuadro 5. Resumen de estadísticos aplicando el FKD a la estación hidrométrica Sayula, considerando la función 
de respuesta modificada. 


Concepto Valor 
Nash-Sutcliffe (Q observado - Q pronosticado) 0.9775 
Nash-Sutcliffe (Q observado - Q actualizado) 0.9808 
Correlación (Q observado - Q pronosticado) 0.9988 
Correlación (Q observado - Q actualizado) 0.9986 
A a) Registros de precipitación efectiva horaria 
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Figura 22. Registros de precipitación y caudal utilizados en el análisis de la cuenca propia de la C.H. Peñitas, incluyendo 
pronóstico de caudales con 1 h de anticipación, utilizando la función de respuesta modificada. 


Caudales “observados” vs. caudales pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-: 2012) 
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Figura 23. Caudales “observados” versus pronosticados en la cuenca propia de la C.H. Peñitas, considerando 
la función de respuesta modificada. 
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Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 24. Caudal “observado” y pronosticado con el FKD (avenida A del 23/11 /2007 01:00 h al 28 /11 /2007 01:00 h), 
aplicado a la cuenca propia de la C.H. Peñitas, considerando la función de respuesta modificada. 
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Figura 25. Caudal observado y pronosticado con el FKD (avenida B del 25/09/2008 03:00 h al 26/09/2008 19:00 h), 
aplicado a la cuenca propia de la C.H Peñitas, considerando la función de respuesta modificada. 
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Caudales “observados” y pronosticados 
Cuenca propia de la C.H. Peñitas. Análisis horario (2005-2012) 
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Figura 26. Caudal observado y pronosticado con FKD (avenida C del 31/10/2009 09:00 h al 03/11/2009 15:00 h), 
aplicado a la cuenca propia de la C.H. Peñitas, considerando la función de respuesta modificada. 


Cuadro 6. Resumen de estadísticos aplicando el FKD a la cuenca propia de la C.H. Peñitas, considerando la función 
de respuesta modificada con pronóstico de 1 h. 


Concepto Valor 
Nash-Sutcliffe (Q “observado” - Q pronosticado) 0.97805 
Nash-Sutcliffe (Q “observado”- Q actualizado) 0.97887 
Correlación (Q “observado”- Q pronosticado) 0.99964396 
Correlación (Q “observado”- Q actualizado) 0.99970262 


de precipitación y 2 h anteriores de caudal, Análisis diario 


calculado con el antitránsito con respecto al 


tiempo k analizado. Los valores de a. y N se 
conservan iguales a los utilizados en los aná- 
lisis realizados anteriormente (figuras 22 a 26 


y cuadro 6). 

La media y desviación estándar de los 
caudales observados fueron de 759.22 y 428.82 
m*/s, respectivamente, mientras que esos 
valores para los gastos pronosticados fueron 
prácticamente los mismos, de 758.97 y 429.61 
m/s. 


Los mejores resultados de este análisis se obtu- 
vieron considerando dos días previos de lluvia 
efectiva y un día de caudal medido dentro de la 
matriz H (figuras 27 a 31 y cuadro 7). 

La media y desviación estándar de los 
caudales observados fueron de 16 972.13 y 10 
246.86 m?/s, respectivamente, mientras que 
esos valores para los gastos pronosticados 
fueron de 16 920.24 y 10 304.58 m*/s. Se 
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Registros de precipitación efectiva horaria 
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Figura 27. Registros de precipitación y caudal utilizados en el análisis de la cuenca propia de la C.H. Peñitas, incluyendo 


pronóstico de caudales con 24 h de anticipación, utilizando la función de respuesta modificada. 


Caudales “observados” vs. caudales actualizados 


Cuenca propia 


Caudales actualizados (FKD) Q (m?/ s) 


de la C.H. Peñitas. Análisis horario (2005-2012) 


Caudales “observados” (calculados con antitránsito) Q (mé 18) 


Figura 28. Caudales “observados” versus actualizados en la cuenca propia de la C.H. Peñitas. Pronóstico con 24 h, considerando 


la función de respuesta modificada. 


preservan los estadísticos aunque en menor 
medida con respecto a los análisis anteriores. 


Discusión de resultados 


El coeficiente de Nash-Sutcliffe para el pro- 
nóstico diario es ligeramente menor que el 


correspondiente al pronóstico horario; sin 
embargo, el intervalo de pronóstico diario es, 
en general, de mayor utilidad que el horario, 
por lo que es aceptable un coeficiente de 
Nash-Sutcliffe menor (aunque aún muy 
alto), a cambio de un mayor intervalo de 
pronóstico. 
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Caudales “observados” y pronosticados 
x104 Cuenca propia de la C.H. Peñitas. Analisis diario (2005-2012) 
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Figura 29. Caudal “observado” y pronosticado con el FKD en la cuenca propia de la C.H. Peñitas (periodo A del 26/9/2006 
al 15/12/2006). Pronóstico con 24 h, considerando la función de respuesta modificada. 
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Figura 30. Caudal “observado” y pronosticado con el FKD en la cuenca propia de la C.H. Peñitas (periodo B del 20/11/2007 
al 09/03/2008). Análisis diario, considerando la función de respuesta modificada. 
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Caudales “observados” y pronosticados 
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Figura 31. Caudal “observado” y pronosticado con el FKD en la cuenca propia de la C.H. Peñitas (periodo C del 27/06/2010 
al 13/01/2011). Análisis diario, considerando la función de respuesta modificada y un pronóstico de 24 h. 


Cuadro 7. Resumen de estadísticos para la aplicación del FKD a la cuenca propia de la C.H. Peñitas, considerando la función 
de respuesta modificada y un pronóstico de 24 h. 
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Concepto Valor 
Nash - Sutcliffe (Q observado - Q pronosticado) 0.83352 
Nash - Sutcliffe (Q observado - Q actualizado) 0.84129 
Correlación (Q observado- Q pronosticado) 0.9996 
Correlación (Q real - Q actualizado) 0.99904 


Coeficiente de Nash-Suttclife 


Coeficiente de Nash-Suttclife 


10.00 


20.00 30.00 40.00 
Intervalos de tiempo en horas 


50.00 


Figura 32. Variación en el Af de pronóstico. 


60.00 


At (h) n nQ Nash-Sutcliffe 
1 12 2 0.97805 
1122 4 2 0.81707 
24 2 1 0.83352 
36 2 2 0.79546 
48 1 1 0.81188 
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Variación del coeficiente de Nash-Sutcliffe en 
términos del intervalo de pronóstico 


En la figura 32 se presentan valores del coefi- 
ciente de Nash-Sutcliffe obtenidos al aplicar 
el algoritmo, considerando distintos At en el 
pronóstico de caudal para la cuenca propia de 
la C.H. Peñitas. 

Es importante señalar que la mayor parte 
de los puntos con mayor diferencia entre los 
valores observados y pronosticados en las 
figuras 23 y 28 se deben a cambios abruptos en 
la señal o bien, en algunos casos, a la ausencia 
de datos. 


Conclusiones 


1. El filtro de Kalman no es un modelo lluvia- 
escurrimiento, sino un algoritmo matemá- 
tico que permite realizar pronósticos del 
estado de un sistema lineal con entradas 
estocásticas. 

2. En el presente trabajo se comprueba que el 
filtro puede ser aplicado de forma exitosa 
en la hidrología, en particular al pronóstico 
de caudales a corto plazo. 

3. La aplicación del filtro de Kalman al pro- 
nóstico de avenidas proporciona resulta- 
dos sumamente aceptables, que permiten 
tener un pronóstico con errores pequeños, 
debido a que éste es corregido con cada 
nueva medición, de manera que podría ser 
considerado una buena opción en el manejo 
y la operación de grandes embalses. 

4. El tiempo computacional invertido es 
mínimo. 

5. Los pronósticos mejorarían de modo consi- 
derable en la medida en que la cuenca esté 
más instrumentada. Así podría evitarse la 
estimación del estado del sistema haciendo 
uso de métodos indirectos, como el anti- 
tránsito de avenidas. 

6. La imposición arbitraria de las variables 
de entrada al algoritmo, como el estado 
inicial y la matriz de la covarianza del error 
inicial, provocan una variación mínima en 
los pronósticos, debido a que el filtro se 


va autocorrigiendo con el paso del tiempo 
hasta disminuir los errores iniciales. 

7. En la mayoría de los puntos en los cuales el 
filtro presenta los errores más significativos 
entre el caudal pronosticado y el observado, 
se presenta un salto abrupto en la señal 
o bien la información de entrada no está 
completa (como cuando alguna estación no 
registra). 

8. Los resultados obtenidos al considerar la 
función de respuesta de la cuenca como 
el hidrograma unitario instantáneo son 
aceptables; sin embargo, la modificación 
propuesta, que incluye tanto los registros 
de caudal como de precipitación en la 
matriz H, mejora de forma considerable los 
pronósticos realizados. 

9. Dentro de las diferentes alternativas utili- 
zadas para definir el valor de la matriz S, el 
mejor resultado se obtiene cuando ésta es 
nula. 

10. Esimportante considerar la posible utilidad 
del FKD en otras áreas (hidráulica, aguas 
subterráneas, etcétera). 
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